Directional Locking Effects and Dynamics for Particles Driven Through a Colloidal 

Lattice 
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We examine the dynamics of a single colloidal particle driven through a colloidal lattice which can 
distort in response to the driven particle. We find a remarkably rich variety of dynamical locking 
phenomena as we vary the angle of the applied drive with respect to the orientation of the colloidal 
lattice. When the driven colloid locks to certain lattice symmetry directions, its motion is not 
necessarily aligned with the drive. Applying a transverse force to the driven particle can result in 
either increased or decreased drag in the driving direction, depending on the angle of the drive. The 
dynamical locking produces anomalies in both the longitudinal and the transverse velocity vs driving 
force curves, including steps and regimes of negative differential resistance. As the interaction of 
the driven particle with the surrounding lattice increases, significant distortion or dislocations in 
the surrounding media occur, and as a result the directional locking is enhanced. We compare these 
results to those obtained for driving particles over fixed substrates, and show that a far richer variety 
of behaviors occurs when the underlying lattice is allowed to distort. We discuss how this system can 
be used for particle species segregation when the onset of different locking angles occurs at different 
drives for varied particle sizes. We also show that the most pronounced locking phases should be 
observable at temperatures up to the melting transition of the colloidal lattice. 

PACS numbers: 82.70.Dd,05.45.-a 



I. INTRODUCTION 

Assemblies of interacting colloidal particles exhibit a 
wide range of ordered, glassy, and liquid states [1,2], and 
have proven an ideal system for studying these equilib- 
rium states with various imaging techniques. Colloids 
can also be studied in nonequilibrium conditions such as 
under shear or when driven with different types of fields 
[2]. Recently Korda et al. [3] investigated the dynamics 
of moving colloids driven over a square two-dimensional 
(2D) optical tweezer array. Here, the optical tweezers 
act as attractive sites [4]. In this system, a series of 
locked transport states occur as the angle at which the 
colloids are driven with respect to the optical lattice is 
varied. In the locked states, the colloids move along a 
high-symmetry direction of the lattice even when this is 
not the same as the direction of drive. The locking pro- 
duces steps in the colloidal velocity vs driving angle, and 
may have applications as a novel method for continuous 
separation of particles of different species, such as macro- 
molecules and biological cells [5-7]. Separation can be 
achieved when a particular species locks to a symmetry 
direction while the remaining species continue to move 
in the direction of the applied drive. Other techniques to 
separate particle species along this line include the trans- 
port of particles through periodic obstacle arrays [8,9]. 

Particle transport through fixed periodic arrays un- 
der varying drive direction has also been investigated in 
simulations of vortices moving in nanostructured super- 
conductors [10] and Josephson junctions [11-13]. Recent 
experiments on vortices driven at different angles through 



periodic pinning arrays find guiding effects similar to the 
features seen in simulations [14]. Similar effects have also 
been studied in electron transport through periodic anti- 
dot arrays [16] as well as atomic friction over ordered sur- 
faces [17]. In general, for all of these systems, the particle 
motion locks to certain symmetry angles of the lattice. 
Refs. [10,16] showed that the velocity vs angle curve has a 
Devil's staircase structure, with a step in the particle ve- 
locity at each locking angle. Certain highly symmetrical 
angles of drive produce stronger locking than others. 

In the experiments of Ref. [3] and the theoretical stud- 
ies of Refs. [10-13,16,17], the periodic lattice through 
which the particles move is completely rigid. Another 
approach which has not been previously considered is for 
the driven particle to move through a stationary, but 
distortable, lattice of particles. In this case, the driven 
particle still interacts with a periodic array of obstacles; 
however, the driven particle can in turn affect the sur- 
rounding media and distort or create dislocations in the 
surrounding lattice. A physical example of such a system 
would be driving a single or small number of colloidal 
particles through a triangular colloidal lattice. This can 
be achieved by placing a colloidal particle that responds 
to an external drive into a colloidal lattice for densities 
and temperatures where the non-driven colloids form a 
triangular lattice. It should be straightforward to drive 
individual colloids with an optical tweezer [7]. It is also 
possible to use magnetic impurity particles driven with a 
magnetic field, as in recent experiments by Weeks et al. 
[18]. A variation on this type of system would be to have 
one colloid or impurity particle fixed at a spot and to 
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drive the colloidal lattice past it. Experiments along this 
line have been performed previously for large obstacles 

A system that has many similarities to colloids is vor- 
tices in superconductors, where the mutual repulsion be- 
tween vortices leads to the formation of a triangular lat- 
tice. For a clean superconductor where the intrinsic pin- 
ning is not strong enough to distort the vortex lattice, it 
may be possible to drive a single vortex with an MFM or 
STM tip. It has also been demonstrated experimentally 
using a lithographically created artificial pin [20] that a 
single pinning site can capture one or more vortices while 
the other vortices are driven past. Other similar systems 
include driving an obstacle through a triangular bubble 
raft lattice, or manipulating a single droplet of a two fer- 
rofluid system in a triangular droplet state. Another sys- 
tem closely related to colloids is 2D dusty plasmas where 
charged dust particles form a triangular lattice [21] and 
individual particles can be driven with lasers. 

In a system where the periodic lattice through which a 
particle is driven can itself be affected by the driven par- 
ticle, a number of intriguing scenarios are possible which 
cannot arise when the substrate is completely rigid. Ex- 
amples include a transition from an ordered flow, where 
the driven particle only distorts the lattice elastically and 
does not create dislocations, to a disordered flow, where 
the particle tears the lattice and dislocations are gener- 
ated. This could occur when the particle moves at an 
angle that is not commensurate with the lattice, or for 
a large enough driven particle. Other effects that could 
occur include the possibility that, for certain driving di- 
rections, the driven particle may drag portions of the 
non-driven lattice with it. Conversely, for other driving 
directions, the particle may be able to slip freely between 
the particles comprising the lattice. Since the lattice can 
locally distort to more easily accommodate the motion 
of the driven particle in the locked regimes, the locking 
effects may be different in nature and more pronounced 
than those observed in previous work on rigid substrates. 
One practical application of this system, as suggested 
in Rcfs. [3,5], is the continuous separation of different 
species of particles as they move through periodic ar- 
rays of obstacles. We demonstrate that the separation 
of particles of different sizes is enhanced by the ability 
of the surrounding lattice to distort. If, under certain 
conditions, the distortable lattice can act as a filter for 
different sized particles, then the additional step of us- 
ing optical tweezer arrays, as in Ref. [3], would not be 
needed. 

II. SIMULATION 

We simulate a 2D monodisperse colloidal assembly us- 
ing Langevin dynamics, with techniques that have been 
described previously [22,23]. The overdamped equation 



of motion for a single colloid i in a system with periodic 
boundaries in the x and y directions is 

r,-^ = t ii +f l d + t r . (1) 

Here 77 = 1 is the damping term arising from the fact 
that the particle is moving through a viscous media. 
The interaction force from the other colloids is fy = 
— SjSi VjV^ry), where the repulsive colloid-colloid in- 
teraction is Yukawa or screened Coulomb, 

V(r ij ) = J ^ n eM-^i-r j \). (2) 

|M T j\ 

Here q^ is the charge of the particle, 1/n is the screen- 
ing length, Yi is the position of particle i, and rj is the 
position of particle j. We initialize the system in a trian- 
gular lattice with TV colloids and a lattice constant a, and 
then insert one additional driven colloid. The screening 
length is fixed at 1/k — 3a. 

The thermal force fr is a randomly fluctuating force 
modeled as random kicks, with the properties < f-r(i) >= 
and < f T (t)f T (t') >= 2r]k B TS(t - t'). For most of the 
results presented here, the colloid density is held fixed 
and the system size is set to L = 34a. With these pa- 
rameters there is a well defined melting temperature T m 
for the unperturbed lattice, which we determine by mea- 
suring the fraction of colloids which are six-fold coordi- 
nated. We measure temperature in units of T m . For the 
first part of the paper we consider the case of T = 
so that the effects of the dynamics can be clearly distin- 
guished. We later show that, for a large range of T/T m , 
essentially the same results appear, and that the lock- 
ing effects are lost for T/T m > 1. Thus the results from 
our work should be applicable to experimental samples 
at finite temperatures. 

We vary the charge on the driven particle qd and the 
direction of drive. All particles except the driven one 
have charge q = 1 and applied dc drive f\ = 0. The force 
on the driven particle has two components: 

f d = .f*St + .f*y. (3) 

We hold f£ fixed at a constant value and increase fjj. 
We then monitor the velocities V x and V y of the driven 
particle. 

We do not take into account possible hydrodynamic ef- 
fects or possible long-range attractions between colloids. 
We have conducted simulations for different system sizes, 
lattice constants a, and screening lengths, and find that 
the qualitative results presented here are robust. 

III. SMALL DRIVEN PARTICLE 

We first consider the case where the driven particle 
interacts only weakly with the surrounding particles so 
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(a) 




(b) " " "(c) 
FIG. 1. Particles (black dots) and trajectories (black lines) 
for a system with q d /q = 0.35 at ft — 1.0 and: (a) 
ft I ft = 0.0, (b) 0.57, and (c) 0.78. 



that, although small distortions of the lattice can occur, 
no dislocations are formed in the surrounding media. We 
define this weak interaction regime as qd/q < 0.5. In the 
case of a triangular lattice aligned in the x-direction, as in 
Fig. 1(a), locking should occur for angles of drive where 
the driven colloid can move freely between the colloids in 
the lattice. This occurs for ratios of the applied drive of 

fy_ = (4) 
fi 2n + l { 1 

where m and n are integers. The corresponding locking 
angle is 9 = arctan(/^ /ft)- Eq. 4 predicts locking at 
60° corresponding to m = 1, n = 0; 30° corresponding 
to m = 1, n = 1; 19.1° corresponding to m = 1, n = 2; 
and 0° for m = 0. Other locking angles occur for higher 
order values of m and n. In Fig. 1 we illustrate the driven 
particle, the surrounding lattice, and the particle trajec- 
tories for a system with fixed parameters of qd/q = 0.35 
and ft = 1.0 for increasing f$/f£ = 0, 0.57, and 0.78. 
The trajectories in all these cases show that the driven 
particle moves through the lattice without significantly 
distorting it, due to the weak charge of the driven par- 
ticle. In Fig. 1(a) we show the commensurate case of 
fy/fx — where the particle channels between the sur- 
rounding colloids. Fig. 1 (b) presents another commensu- 
rate case for fy/ft = 0.57 where the particle moves in a 
periodic orbit along a symmetry line of the lattice. This 
orbit corresponds to m — 1, n = 1 from Eq. 4, so the 
flow is along 30°. In Fig. 1(c) we show an incommensu- 
rate case of fy/f^ = 0.78, which is not a simple 
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FIG. 2. The power spectra obtained from a time series of 
V y for: (a) the commensurate case fy/fx = 0.57 shown in 
Fig. 1(b), and (b) the incommensurate case fy/fx = 0.78 
shown in Fig. 1(c). 

ratio in Eq. 4. Here the particle does not move in a 
completely periodic orbit but shows some non-periodic 
deviations; however, the particles in the surrounding lat- 
tice still show little deviation from their positions. We 
observe similar types of motion at the other incommen- 
surate angles as well. 

In order to better demonstrate the difference between 
the commensurate and incommensurate particle orbits in 
Fig. 1, we analyze the power spectra of the velocity V y 
of the driven particle for the two cases. In Fig. 2 we 
show the power spectra S(y) obtained from a time series 
of V y of the driven particle at fixed fy/ft- Figure 2(a) 
illustrates S(y) for the commensurate case fy/ ft = 0.57 
shown in Fig. 1(b). Here the spectrum has a prominent 
peak at v = 0.0127 inverse molecular dynamics steps with 
higher harmonics, indicating that the particle motion is 
strictly periodic with a single frequency. In a commen- 
surate orbit, a single frequency is expected to appear 
since the particle is moving along a symmetry direction 
of the lattice and is slowed by the interactions from the 
surrounding lattice at a constant spacing and a constant 
rate. In Fig. 2(b), the spectrum for the incommensu- 
rate case fy/ft = 0.78 shown in Fig. 1(c) has many 
peaks, indicating that the particle is undergoing motion 
with many different periods giving a broad spectrum. 
In general, for these incommensurate phases, the par- 
ticle likely jumps between several different closely spaced 
locked phases with high m and n values which, together, 
produce an angle of motion which is close to the angle of 
the drive. We find similar behaviors in the power spectra 
at other commensurate and incommensurate angles. 

In Fig. 3(a) we simultaneously plot V x and V y vs ft /ft, 
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FIG. 3. (a) V x (dark line) and V y (light line) for fixed 
= 1.0 and increasing for q d /g = 0.35. (b) V y (dashed 
line) and its derivative R y (solid line), showing the locking 
phases more clearly. 

with fixed f~ — 1.0, showing that the locked phases 
strongly affect the transport proprieties. We highlight 
a few ratios of (to, n) from Eq. 4 to indicate the most 
prominent locking phases: (0, 0), which starts at fy/f x = 
in the bottom left corner of Fig. 3(a); (1, 1) for flow at 
30°; and the very prominent (1,0) locking centered at 
fy/fx — T7. There are additional smaller steps that 
are difficult to see for this value of q^. Another inter- 
esting feature is that V x has an average value of 0.9 for 
fy/fx < 2-0. We note that, for a single driven colloid 
moving in the absence of any other particles, we would 
have V x = 1.0, which is higher than the value obtained 
in the presence of the other particles. This indicates that 
the driven particle experiences an additional damping 
due to the interactions with the other particles. This 
damping originates when the surrounding colloids shift 
slightly as the driven particle moves past. Since the mo- 



tion of these surrounding colloids is also overdamped, the 
only source of energy is the driven particle. In the case 
of a rigid substrate [10], this additional damping effect 
is absent since the periodic substrate cannot distort and 
absorb energy. Another immediately obvious difference 
between this system and particle motion over a rigid sub- 
strate is that in Fig. 3(a), none of the steps in V y have 
dV y /df£ = [1,10,16]. Although the slope of V y is re- 
duced on the steps, it does not drop to zero, indicating 
that along the locking directions the particle velocity V y 
continues to increase with f£. Fig. 3(a) also indicates 
that, in the locking regions, V x is not constant but shows 
features where the velocity both increases and decreases 
which coincide with slope changes in V y . On the steps, 
when V y locks to the lower slope value, V x locks to a 
constant positive slope which can be seen most clearly 
near f*/f* = 1.7 at the 60° locking of (1,0). At the 
lower edge of the (1, 1) and (1,0) steps, V x decreases in 
an abrupt dip before rising linearly on the step. On the 
upper end of these steps, when the slope in V y increases 
again, V x shows another decrease. 

In Fig. 3(b) we plot V y along with its derivative R y to 
show the locking phases more clearly. At the (1,0) lock- 
ing, R y drops to a lower value due to the smaller slope 
of V y in the locking regions. Peaks appear in R y at the 
beginning and end of each locked phase, reflecting the in- 
creased slope of V y outside of the locked regions. R y also 
shows some oscillations above the (1,0) locking, which 
we discuss in more detail in Section IV-C. Additionally, 
locking phases appear as dips in R y for the (1,1) locking 
centered at fy/f x = 0.57 as well as at (1, 2) and (1, 3). 

Fig. 3(a) also demonstrates that there is a clear trans- 
verse depinning threshold for the (0, 0) locking: V y = 0.0 
for fy/f x < 0.15. Here, despite the fact that the particle 
is moving in the a;-direction, there is an effective pinning 
threshold for motion in the y-direction. Below this de- 
pinning threshold, the particle trajectories are the same 
as those shown in Fig. 1(a), where the particle moves in 
an effective ID trough. Transverse depinning thresholds 
have been observed in various systems of particles driven 
longitudinally over a 2D periodic substrate as an increas- 
ing transverse force is applied. These systems include 
longitudinally moving vortices in periodic pinning arrays 
[10,11,13] and atoms moving over periodic surfaces [17]. 
In colloidal experiments on periodic substrates, a locking 
is observed for small angles [3]. In all these cases the 
particle motion remains locked in the longitudinal direc- 
tion until a large enough transverse force is applied. In 
Fig. 3(a) as fy/f x approaches the transverse depinning 
threshold, V x , the particle velocity in the ^-direction, 
decreases. This decrease occurs because the increasing 
transverse force shifts the driven particle closer to the 
lattice particles, increasing the interaction between the 
driven particle and the other particles. Since the sur- 
rounding lattice is flexible, some of the motion of the 
driven particle is transfered into small 



4 




1 2 



FIG. 4. V x vs f*/f* for q d /q = 0.125, 0.25, 0.35, and 0.5 
(from top left to bottom left). 

distortions of the surrounding lattice. These distortions 
increase as the driven particle is forced closer to the par- 
ticles in the lattice. Once the driven particle begins to 
move in the y-direction, its velocity in the x-direction 
increases since the driven particle no longer approaches 
the particles in the surrounding lattice as closely and thus 
experiences less damping of its motion. Remarkably, for 
1.8 < f*/f£ < 2.0, V x is greater than 1.0, indicating that 
the particle is moving faster than it should given the ap- 
plied driving force in the x direction. This indicates that 
some of the energy from the y component of the drive 
is being transferred to the x-direction. Above the (1,0) 
step for fy/fx > 2.0, V x decreases slowly back to 1.0. 

We next consider how the locking regions evolve with 
increasing q d . In Fig. 4 we plot V x vs f$/f£ for q d /q = 
0.125, 0.25, 0.35, and 0.5 with f* = 1.0. The plots are 
overlaid with no offset. For fy/f£ — 0, V x decreases for 
increasing qd, indicating that the damping increases as qd 
rises. This decrease in V x occurs because the driven par- 
ticle interacts more strongly with the surrounding parti- 
cles, causing larger distortions and transferring more of 
its energy into the lattice. Additionally, all of the lock- 
ing phases become wider with increasing qd, indicating 
that as the distortions in the surrounding lattice from the 
driven particle increase, the channeling effect becomes 
stronger. The locking phases for the higher order values 
of (m, n) can now be clearly resolved, especially for the 
qd/q = 0.5 curve where a particularly large number of 
dips occur for fy/f x < 1.0. For the locking at (1, 0), the 
value of V x exceeds 1.0 for fy/f x > 1.7 at all values of qd- 
Just above the (1, 0) locking, V x decays back to 1.0. Dur- 
ing this decay, there are clear small scale periodic peaks 
which are a real effect that we study later in this paper. 

We note that for small qd, the behavior of this system 
is similar to that of a particle moving through a rigid 




FIG. 5. (a) V y vs fy/fx at the transverse depinning tran- 
sition for the system shown in Fig. 3 with qd/q — 0.35. (b) 
Log-log plot of V y vs (fy/fx - fc) where f c = 0.15 is the 
transverse depinning threshold. The solid line is a power law 
fit with an exponent of /3 = 0.58. 



lattice in Ref. [10], where the same type of anomalous 
steps were observed in both V x and V y . The locations 
of the steps in Ref. [10] occurred for different values of 
fy/fx since the underlying lattice was square rather than 
triangular. Additionally, the slope of V y along the steps 
was nearly zero in the rigid lattice, making V y appear 
much more step-like. The similar behavior found here at 
small qd occurs when the driven particle can no longer 
cause significant distortions in the surrounding lattice, 
rendering it effectively rigid. 

A. Dynamics of Transverse Depinning 

We now examine the initial transverse depinning tran- 
sition in more detail. We first consider the shape of V y vs 
fy/fx near the transverse depinning. For elastic media 
depinning in the presence of quenched disorder, theory 
predicts V oc (f d — f c )^ , where f c is the depinning thresh- 
old [24] . For a single particle depinning from a sinusoidal 
substrate in ID, (3 — 1/2. In 2D systems of collectively 
interacting particles where plastic deformation occurs at 
depinning, studies find scaling in the velocity force curves 
with f3 > 1.5 [25]. It is not known if there would be any 
scaling for the transverse depinning of a longitudinally 
moving particle. In Fig. 5(a) we show V y near the onset 
of motion in the y-direction for the system in Fig. 3 with 
q d /q = 0.35. We plot V y vs fy/f x - f c on a log-log scale 
in Fig. 5(b), where f c = 0.15 is the transverse depinning 
force. We find a reasonable power-law fit with j3 = 0.58 
(solid line in Fig. 5(b)), close to the single particle value 
of j3 = 1/2. This indicates that the transverse depinning 
is elastic in nature, and that when the particle depins it 
does not induce tearing in the surrounding lattice. 

For the steps above (0, 0), whenever the particle leaves 
a locked region, V y shows curvature similar to that of the 
initial transverse depinning. This suggests that on the 
higher order steps, when the particle is channeling along 
an easy flow direction, the particle is effectively pinned in 
the direction transverse to td, rather than in the direction 
transverse to f£. As f$ increases, the transverse force 
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(c) (d) 

FIG. 6. Particles (black dots) and trajectories (black lines) 
near the transverse depinning transition for the system shown 
in Fig. 5 at f*/f c = (a) 1.05, (b) 1.22, (c) 1.34, and (d) 1.5. 

eventually becomes high enough for the particle to depin 
and begin moving at a new angle. 

In Fig. 6 we plot the particle trajectories at the trans- 
verse depinning transition for the system shown in Fig. 5. 
In Fig. 6(a), just above the initial transverse depinning 
at fy/fc = 1-05, the particle stays locked along the x- 
direction for about 16a before moving in the positive y 
direction to the next row. The particle moves in a stair- 
case fashion. In Fig. 6(b) for fy/fc = 1-22, we find a 
staircase motion similar to that of Fig. 6(a), but with 
the particle translating in the y direction every 7a, pro- 
ducing a net V y that is higher than that of Fig. 6(a). In 
Fig. 6(c) at f*/f c = 1.34 and Fig. 6(d) at f d Jf c = 1.5, 
the staircase motion persists with the particle moving 
shorter distances in the x-direction between jumps to the 
next row in the y direction. 



IV. INTERMEDIATE SIZED DRIVEN 
PARTICLES 

Next we consider the case for 0.5 < qd/q < 3.0, which 
we term the intermediate particle regime. Here we ob- 
serve that the locking phases become more pronounced 
and appear more step- like in V y with dV y /dfy — in some 
places. We also find regimes where V y decreases for in- 
creasing fy, indicating a negative differential transverse 
resistance, or dV y /dfy < 0. The scaling in the initial 
transverse depinning is also lost and is replaced by 
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FIG. 7. (a) V x (dark line) and V y (lighter line) vs f$/f* 
for q d /q = 1.0. (b) Closeup of (a) for 0.0 < f$/f£ < 1.0 

a sharp jump, indicating that for qd/q > 0.5 the system 
is no longer in the elastic regime. In Fig. 7(a) we show V x 
and V v for a system with qd/q = 1.0. For these interme- 
diate qd values, the width of the (1,0) locking becomes 
so large that it begins to overlap with the smaller lock- 
ing phases, and some of the smaller locking steps are lost. 
The dips in V x are more pronounced and there are several 
regions where both V x and V y show jumps into and out of 
different locked phases, rather than the smoother transi- 
tions seen for lower values of qd- The slopes of V x and V y 
along the (1, 0) locking are both positive but show a non- 
linear bowing effect. In Fig. 7(b) we plot V x and V y for 
< fy/ f x < 1, showing the smaller locking phases more 
clearly. Along the step (1,1), centered at f y /f x = 0.6, 
the slope of V y shows a bow feature, indicating that the 
slope becomes negative near fy/f x = 0.73. On the next 
step we find a similar bow feature. This behavior is in 
contrast to the studies with rigid substrates, where the 
slope of V y never drops below zero along the 
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FIG. 8. Particles (black dots) and trajectories (black lines) 
for the system in Fig. 7 with f$/f£ = (a) 0.39, (b) 0.6, (c) 
0.8, and (d) 2.0. 



steps. Fig. 7 also shows that both V x and V y can de- 
crease simultaneously, indicating an overall increase in 
the damping. For these intermediate qd values there is 
considerably more distortion in the surrounding lattice. 
Just before the driven particle exits a locked region, it 
has the strongest interaction with the surrounding me- 
dia and experiences the largest damping. At the jumps 
in velocity accompanying each step, the driven particle 
can cause enough distortion to temporarily generate a 
localized dislocation in the lattice. 

A. Trajectories in the Locked Phases 

In Fig. 8 we illustrate some of the particle trajectories 
for different locking regimes from the system in Fig. 7. 
Fig. 8(a) shows the step near fy/f£ — 0.39, which corre- 
sponds to the m = 1, n = 2 locking phase. The particle 
moves in a zig-zag pattern. The elastic distortions in 
the surrounding colloid lattice arc clearly visible in the 
form of very small loops on both sides of the path of the 
driven particle. Here, the colloids closest to the driven 
particle move a distance less than a as the driven particle 
approaches, and then return to their initial positions af- 
ter the driven particle passes. There are some additional 
smaller distortions in the lattice at larger distances from 
the driven colloid; however, these motions are too small 
to be visible at the resolution in Fig. 8. In general, when 
the intermediate sized particles move through a channel 
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FIG. 9. V y vs fy/fx for varied qa/q- A systematic shift 



is added in the y direction for clarity. From top to bottom, 
q d /q = 2.0, 1.75, 1.5, 1.0, and 0.75. 

during a locked phase, the motion is elastic and no dislo- 
cations are generated. Only during the transitions out of 
the locked phases do temporary distortions greater than 
a in the surrounding lattice occur. 

Fig. 8(b) shows the trajectories at fy/f£ = 0.6, when 
the particle moves at 30° in the (1, 1) locking regime. The 
width of this (1,1) locked region is considerably larger 
than that of the (1, 2) phase in Fig. 8(a). Here, the par- 
ticle follows a sinusoidal path and causes distortions in 
the surrounding lattice; however, these distortions are of 
a smaller magnitude than those in Fig. 8(a), indicating 
that the particle channels more easily on the (1, 1) step. 
In Fig. 8(c) we plot the trajectories for the step region 
at fy/fx = 0.8, where more significant distortions in the 
surrounding lattice occur than for Fig. 8(b). The width 
of this step is smaller than the (1, 1) step in Fig. 8(b). In 
Fig. 8(d) we show the trajectories at fy/f£ = 2.0 for the 
most prominent step, (1,0), where the vortices channel 
at 60°. There are almost no distortions in the surround- 
ing lattice. For a value of fy/f^ just below the end of 
the (1,0) phase, fy/f£ ~ 2.6, significant distortions of 
the colloids to the upper left of the driven particle oc- 
cur in the surrounding lattice, since the driven particle 
is pushed in this direction by the increasing fy. When 
these distortions become large enough, the (1,0) phase 
ends. 

B. Locked Phase Evolution as a Function of qd 

In Fig. 9 we plot V y vs fy/ fx for increasing values of qd- 
From bottom to top, qd/q = 0.75, 1, 1.5, 1.75, and 2.0. In 
this range of qd, the initial transverse depinning threshold 
changes little; however, the width of the (1, 0) locking 
which begins at ft/ ft ~ 1.1 grows until it saturates 



7 



fV 

y x 

FIG. 10. 14 vs f*/fg for the system in Fig. 9 with a sys- 
tematic shift added to the y direction for clarity. From top to 
bottom, q d /q = 2.0, 1.75, 1.5, 1.0, and 0.75. 

at qd/q = 1.75. The width of the locking regions for 
fy/fx < 1-1 do not increase much with qd; however, 
they became increasingly smeared, which is particularly 
noticeable for qd/q > 1.75. The smearing occurs when 
the lattice distortions become large enough to create dis- 
locations in the surrounding lattice. These dislocations 
appear when the particle jumps from one symmetry an- 
gle to another, and thus they smear out the transition. 
The saturation of the width of the main (1, 0) step also 
coincides with the onset of the creation of dislocations as 
the particle passes into and out of the locked phase. In 
Fig. 10 we plot the corresponding V x vs fy/fx for the 
system in Fig. 9. Again we find a smearing of the fea- 
tures with increasing qd for fy/fx < 1-0. Additionally, 
we find an increase in the magnitude of the fluctuations 
in V x just above the end of the (1, 0) locking phase for 
qd/q = 2.0, due to the formation of dislocations in the 
surrounding media. 

In Fig. 11 we plot a phase diagram showing the evolu- 
tion of the widths of the transverse pinned phase T p and 
the (1,0) locked phase L. The phases are identified from 
the features in V y shown in Fig. 9. Here the transverse 
pinned phase T p grows with qd until its width saturates 
around qd/q = 0.75, while the locked phase L also grows 
with qd and saturates at qd/q = 1-75. There is a small de- 
crease in the width of the locked phase for higher values 
of qd. Above qd/q — 2.5, the phases become increas- 
ingly difficult to distinguish due to the large distortions 
in the lattice that cause large velocity fluctuations near 
the transitions into and out of the locked phases. 

The (1,0) phase grows with increasing particle inter- 
action strength qd because the driven particle can more 
easily distort the surrounding lattice, allowing the lat- 
tice to accommodate the (1,0) phase for larger angles 
between the driving force and 60°. Eventually, for large 



d d 

y x 

FIG. 11. Phase diagram for qd/q vs fy/fx- The transverse 
pinned phase is marked T p , and the locked phase for m = 1, 
n = is marked L. 

enough qd/q, the distortion in the surrounding lattice 
become strong enough to allow dislocations to be created, 
enabling the particle to jump out of the locked channel, 
and limiting the growth of the locked phase L in Fig. 11. 
As the particle interaction increases above qd/q — 2, this 
transition occurs at lower values of fy/fx, and the width 
of the (1,0) phase decreases as seen in Fig. 11. 

The behavior of the locked phase in Fig. 11 also indi- 
cates how this system could be used for particle species 
separation. If two different species of particles are moving 
through the sample with the same applied fy/fx, then 
it is possible for one of the species to be in the locked 
L phase while the other species is not, and thus the two 
species will move in different directions. For example, for 
a drive of fy/fx = 1-25, particles with qd/q > 0.75 are 
in the L phase while particles with qd/q < 0.75 are not 
and will move at smaller angles. One advantage of this 
technique over the optical traps is that particles with the 
same optical properties but different charges or sizes 




FIG. 12. (a) V x - Vo vs fi/f£ - f c , where Vo = 1.0 and 
f c = 2.338 is the value of the driving force when the (1,0) 
locking phase ends for a system with qd/q = 0.75. (b) 
Log-linear plot of (a). The solid straight line is a fit with 
the form Aexp(-bx) where b = 2.0 and A — 0.1. 




FIG. 13. (a) V y vs for the system in Fig. 12, just 

after the (1,0) locking ends, (b) Corresponding V x . 



could be separated. 

C. Dynamics just outside of the m — 1, n — phase 

We now look more closely at some features in the V x 
and V y curves as well as the particle dynamics just at the 
end of the m = 1, n = phase. In Fig. 10, after the end 
of the (1, 0) phase, the V x curves decay back to V x = 1.0. 
All the values qd/q < 2.5 show a similar decay in this 
regime. To examine this decay, in Fig. 12(a) we plot 
V x - V vs fg/f* - f c , where q d /q = 0.75, V = 1.0, and 
f c = 2.338 is the driving force at which the (1, 0) locked 
regime ends and the V x decay begins. In Fig. 12(b) we 
plot the same curve in log-linear form and show a good 
fit to an exponential form Aexp(-bx), with b = 2.0. For 
other values of qd/q we observe the same behavior with 
similar values for 6. The fitting constants and curve shape 
do not change when the drive is swept more slowly in 
the simulations, indicating that we are not observing a 
transient behavior. 

In the decay region shown in Fig. 12, there appear to 
be some smaller oscillations in V x . Similar oscillations 
are also seen in V y . Focusing on this region, we find that 




FIG. 14. Particles (black dots) and trajectories (black 
lines) for the system in Fig. 13 at: (a) the step at 
ft/ ft = 2.412, and (b) the step at = 2.51. 



the oscillations are actually steps. In Fig. 13(a) we plot 
V y for the case of q d /q = 0.75 right after the end of the 
(1,0) locked phase. Here V y does not increase smoothly 
but passes through a series of small jumps. Between these 
jumps V v increases linearly. In Fig. 13(b) we plot the cor- 
responding V x for the same range of driving forces. Here 
V x also decreases in a series of jumps which coincide with 
the jumps in V y . Between the jumps, V x increases lin- 
early. For higher values of fy/f x beyond what is shown 
in Fig. 13, the jumps continue and become more frequent 
until they overlap. In order to determine the cause of the 
steps, in Fig. 14 we plot the particle trajectories on the 
small step centered at fy/f x = 2.412 from Fig. 13. Here 
the particle moves at 60° for a distance of about 17a, 
then jumps over to the next row, moves at 60° again, and 
then repeats the cycle of alternately jumping and mov- 
ing. In Fig. 14(b), on the step just after f*/f* = 2.51 
a similar motion is seen as in Fig. 14(a); however, the 
particle moves a distance of about 9a before jumping to 
the adjacent row. This staircase type motion is similar to 
the dynamics just above the transverse depinning transi- 
tion as shown in Fig. 6. In both these cases the particle 
moves predominantly along the symmetry angle which is 
0° for the transverse depinning and 60° here, with peri- 
odic jumps into the next row. For the higher order steps, 
the particle moves a shorter distance in the 60° direction 
before jumping to the next row. 

V. STRONGLY INTERACTING PARTICLES 

We next consider the case of driven particles with 
q<i/q > 2.5. In this regime we find that a significant 
number of dislocations are introduced into the system. 
In Fig. 15 we show V y and V x for qd/q — 5.0. Here 
the curves are noisy with several sharp steps. No clearly 
defined locking regime is present; however, there are rem- 
nants of the (1,0) locking regime. 

The smaller phase locking regions for drive angles less 
than 60° that were present at smaller qd are completely 
washed out in Fig. 15. On these smallest steps, as 
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FIG. 15. (a) V y vs for q d /q = 5. (b) Corresponding 

•d j fd 
y / J x ■ 



V x vs f$/f? 



illustrated in Fig. 8, the driven particle created the 
largest distortions in the surrounding lattice. As qa is 
increased, dislocations appear first for the steps which 
produce the largest lattice distortions, and hence the 
smallest phase locking regions are smeared out first. This 
explains the disappearance of these phases at large qd- 
Fig. 15 shows that there is still an appreciable locking 
effect for the m = 1, n = phase. The intermittent 
jumps in the velocities are due to the formation of dis- 
locations, a process which abruptly changes the velocity 
of the driven particle. Additionally, a small number of 
dislocations remain in the colloid lattice for fy/f£ below 
the m = 1, n = phase. In Fig. 16 we show the trajecto- 
ries for the system in Fig. 15 for fy/f£ = 1.0. Here the 
trajectories are strongly disordered and in general do not 
overlap. For this value of drive a number of dislocations 
are created in the lattice. Some of the dislocations dis- 
appear when the lattice reorders once the driven particle 
moves past; however, some persist for long times. Since 
we are using periodic boundary conditions, the particle 
passes through the same region repeatedly and interacts 
with the persistent dislocations, which could make the 



FIG. 16. Particles (black dots) and trajectories (black 
lines) for the system in Fig. 15 at fy/ft = 10. 

flow more disordered. To address this we have performed 
the same simulations for larger systems and examined the 
trajectories before the particle comes back on its previous 
path. We find the same disordered flow, indicating that 
it is the process of dislocation creation which causes the 
particle trajectories to become disordered, and not the 
existence of dislocations from a previous passage. We 
have also measured the power spectra in this case and 
find white noise rather than a narrow band signal. 



VI. TEMPERATURE EFFECTS 

In a real colloidal system, temperature effects will be 
relevant. For a triangular colloidal system at a finite 
density, there is a well defined melting temperature T rn 
above which dislocations proliferate. We investigate the 
effects of the temperature on the locking up to T m . We 
concentrate on the system with qd/q = 0.5 where locking 
is still clearly visible but the driven particle does not 
generate dislocations in the surrounding media at T = 0. 
In Fig. 17(a) we show V x for T/T m = 0.0, 0.44, 0.7, 0.9, 
and 1.01. For the lowest temperature, the width of the 
(1,0) step is the largest. It decreases for higher T and 
disappears above T m . Additionally, the initial dip in V x 
for the transversely pinned phases is also reduced as the 
temperature increases. The higher order locked regions 
for fy/fx < 1A become washed out for T/T m > 0.7. 

In Fig. 17(b) we plot the width W of the (1,0) lock- 
ing step from Fig. 17(a). Here the width deceases with 
increasing T and drops suddenly at T m . This indicates 
that, even for relatively small particles, the locking effects 
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FIG. 17. (a) V x vs f$/f£ for varied T/T m = 0.0, 0.44, 0.7, 
0.9, and 1.01. T m is the temperature at which the non-driven 
lattice melts. The curve with the largest (1,0) locking phase 
(centered at f$/f£ = 1.7) is T = 0.0. The width of the (1, 0) 
locking phase decreases for increasing T. (b) Width W of the 
(1,0) locking vs T/T m for the system in Fig. 17(a). 

should be observable in an experimental system all the 
way up to the melting transition of the colloidal lattice 
itself. 



VII. SUMMARY 

In summary, we have investigated the dynamics of 
driven particles moving through a triangular colloidal 
lattice for varied orientations of the drive with respect 
to the symmetry of the colloidal lattice. This system, 
where the driven particle causes distortions in the sur- 
rounding colloidal lattice, differs from previously studied 
systems in which particles are driven through a lattice 
of fixed obstacles or traps. We find a series of locked 



phases as a function of drive angle where the particles 
prefer to move along certain symmetry directions. For 
driven particles that interact only weakly with the sur- 
rounding lattice, we find a Devil's staircase structure. For 
larger driven particles, the locking effects become much 
more pronounced and the transitions into and out of the 
locked states become much sharper. These transitions 
can be accompanied by the creation of short-lived dis- 
locations in the surrounding lattice. We also find new 
features such as negative differential resistance in both 
the transverse and longitudinal velocities in the interme- 
diate particle regime. As the driven particles increase in 
size, the smallest locking regimes are lost when the main 
locking regimes grow and significant distortion of the sur- 
rounding lattice occurs. For larger particles, the locking 
effects are washed out when significant numbers of dislo- 
cations are created. Since the width of the locking steps 
depends strongly on the size of the driven particle, our 
results may be useful as a method for particle separation 
techniques or electrophoresis. We have also studied the 
dynamics at the transitions into and out of the main lock- 
ing states and find a staircase-like motion of the driven 
particle. 
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useful discussions. This work was supported by the U.S. 
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